Preference of Small Molecules for Local Minimum Conformations when Binding to Proteins

It is well known that small molecules (ligands) do not necessarily adopt their lowest potential energy conformations when binding to proteins. Analyses of protein-bound ligand crystal structures have reportedly shown that many of them do not even adopt the conformations at local minima of their potential energy surfaces (local minimum conformations). The results of these analyses raise a concern regarding the validity of virtual screening methods that use ligands in local minimum conformations. Here we report a normal-mode-analysis (NMA) study of 100 crystal structures of protein-bound ligands. Our data show that the energy minimization of a ligand alone does not automatically stop at a local minimum conformation if the minimum of the potential energy surface is shallow, thus leading to the folding of the ligand. Furthermore, our data show that all 100 ligand conformations in their protein-bound ligand crystal structures are nearly identical to their local minimum conformations obtained from NMA-monitored energy minimization, suggesting that ligands prefer to adopt local minimum conformations when binding to proteins. These results both support virtual screening methods that use ligands in local minimum conformations and caution about possible adverse effect of excessive energy minimization when generating a database of ligand conformations for virtual screening

Normal-Mode-Analysis–Monitored Energy Minimization Procedure for Generating Small–Molecule Bound Conformations

The energy minimization of a small molecule alone does not automatically stop at a local minimum of the potential energy surface of the molecule if the minimum is shallow, thus leading to folding of the molecule and consequently hampering the generation of the bound conformation of a guest in the absence of its host. This questions the practicality of virtual screening methods that use conformations at local minima of their potential energy surfaces (local minimum conformations) as potential bound conformations. Here we report a normal-mode-analysis–monitored energy minimization (NEM) procedure that generates local minimum conformations as potential bound conformations. Of 22 selected guest–host complex crystal structures with guest structures possessing up to four rotatable bonds, all complexes were reproduced, with guest mass–weighted root mean square deviations of <1.0 Å, through docking with the NEM–generated guest local minimum conformations. An analysis of the potential energies of these local minimum conformations showed that 22 (100%), 18 (82%), 16 (73%), and 12 (55%) of the 22 guest bound conformations in the crystal structures had conformational strain energies of less than or equal to 3.8, 2.0, 0.6, and 0.0 kcal/mol, respectively. These results suggest that (1) the NEM procedure can generate small–molecule bound conformations, and (2) guests adopt low-strain–energy conformations for complexation, thus supporting the virtual screening methods that use local minimum conformations

Peptide Bond Distortions from Planarity: New Insights from Quantum Mechanical Calculations and Peptide/Protein Crystal Structures

By combining quantum-mechanical analysis and statistical survey of peptide/protein structure databases we here report a thorough investigation of the conformational dependence of the geometry of peptide bond, the basic element of protein structures. Different peptide model systems have been studied by an integrated quantum mechanical approach, employing DFT, MP2 and CCSD(T) calculations, both in aqueous solution and in the gas phase. Also in absence of inter-residue interactions, small distortions from the planarity are more a rule than an exception, and they are mainly determined by the backbone ψ dihedral angle. These indications are fully corroborated by a statistical survey of accurate protein/peptide structures. Orbital analysis shows that orbital interactions between the σ system of Cα substituents and the π system of the amide bond are crucial for the modulation of peptide bond distortions. Our study thus indicates that, although long-range inter-molecular interactions can obviously affect the peptide planarity, their influence is statistically averaged. Therefore, the variability of peptide bond geometry in proteins is remarkably reproduced by extremely simplified systems since local factors are the main driving force of these observed trends. The implications of the present findings for protein structure determination, validation and prediction are also discussed